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THE STABILITY OF A TRAILIN G-LINE VORTEX 
IN COMPRESSIBLE FLOW 


Jillian A. K. Stott and Peter W. Duck 1 

Department of Mathematics 
University of Manchester 
Manchester, M13 9PL 
UNITED KINGDOM 

ABSTRACT 

We consider the inviscid stability of the Batchelor (1964) vortex in a compressible flow. 
The problem is tackled numerically and also asymptotically, in the limit of large (azimuthal 
and streamwise) wavenumbers, together with large Mach numbers. The nature of the solution 
passes through different regimes as the Mach number increases, relative to the wavenum- 
bers. At very high wavenumbers and Mach numbers, the mode which is present in the 
incompressible case ceases to be unstable, whilst a new “centre mode” forms, whose stabil- 
ity characteristics are determined primarily by conditions close to the vortex axis. We find 
that generally the flow becomes less unstable as the Mach number increases, and that the 
regime of instability appears generally confined to disturbances in a direction counter to the 
direction of the rotation of the swirl of the vortex. 

Throughout the paper comparison is made between our numerical results and results 
obtained from the various asymptotic theories. 


! This research was supported by the National Aeronautics and Space Administration under NASA Con- 
tract Nos. NAS1-18605 and NAS1-19480 while the second author was in residence at the Institute for 
Computer Applications in Science and Engineering (ICASE), NASA Langley Research Center, Hampton, 


VA 23681. 
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1. Introduction 


In recent years there has been a good deal of interest in the stability of incompressible 
swirling vortex-type flows. Two important applications to this area of research are the 
breakdown of trailing-line vortices behind aircraft and to tornadoes; this class of flow may 
also be applicable to flows inside turbines and compressors, to which the present work 
would be particularly relevant. 

The earliest works in the area of the stability of swirling vortex flows include those 
of Lessen & Paillet (1974) and Lessen, Singh & Paillet (1974). In the former paper the 
stability of the Batchelor (1964) vortex was considered, at finite Reynolds numbers up to 
150. In the second paper, the inviscid stability of this vortex was studied and revealed an 
increase in growth rate at increasingly large wavenumbers (with disturbances being most 
dangerous counter to the direction of the swirl). 

Duck & Foster (1980) showed that for a given wavenumber a multiplicity of modes 
exists. Leibovich & Stewartson (1983) and Duck (1986) considered the limit of large 
wavenumber for this problem, and showed that a finite (maximum) growth rate was at- 
tained. The aforementioned studies suggested an upper and lower neutral value of axial 
wavenumber. The upper neutral point for large azimuthal wavenumber was treated by 
Stewartson & Capell (1985), who showed that the “ring mode” structure of the unstable 
modes persisted near the upper neutral points. Stewartson & Brown (1985) considered 
these upper neutral points for order one azimuthal wavenumbers, and found that in this 
case the modes were of centre mode type, similar to those found in a related study on 
swirling Poiseuille flow (Stewartson & Brown 1984). The behaviour of the unstable modes, 
close to the lower neutral point, at large azimuthal wavenumber, was investigated by Stew- 
artson & Leibovich (1987), and determined that in this case the instability disturbances 
were centred near the axis of the vortex. 

More recently, viscous results have been presented at finite (but large) Reynolds 
numbers by Khorrami ei al. (1989) and Khorrami (1991). In this latter paper, it was 
shown that additional unstable modes exist, in which viscosity plays a destabilising role. 
These modes were analysed by Duck & Khorrami (1991). The inviscid analysis is also 
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applicable to other vortex profiles, including that of Long (1961), as studied by Foster & 
Duck (1982) and Foster Sz Smith (1989). 

Little attention has been paid to the stability of compressible vortex flows of this 
type. On the other hand the area of compressible jet flow has been investigated for some 
years now, the work of Michalke (1971, 1984) being relevant here, although restricted 
to non-swirling flows. Compressible swirling jet flows have also received some attention, 
the work of Coleman (1989) should be mentioned, who studied the superposition of a 
Rankine vortex on a top-hat jet velocity field. More recently, Khorrami (1991) studied 
a compressible swirling axisymmetric jet, by assuming the incompressible flow of Gortler 
(1954) and Loitsyanskii (1953) was applicable in the compressible regime. 

In this paper, we take cylindrical polar coordinates, ( Ir , 9, lx), with the x axis lying 
along the axis of the vortex (which is taken to be axisymmetric), and l is some stream- 
wise scale. We also take the flow far from the vortex centre to be directed along the 
x-direction. The velocity field is written as U*u = U*(u, v, w), the fluid density is p^p, 
temperature T^T, first and second coefficients of viscosity p^p, pZo^ respectively, and 
pressure p^U^p. Here superscript asterisk and subscript oo denote dimensional and 
freestream variables, respectively, and also U* is a velocity scale, defined below in ( 2 . 6 ), 
whilst UZc is the freestream velocity. We define the flow Reynolds number 


Re = 


PtoW 


( 1 . 1 ) 


and we have a flow Mach number given by 
7T* 

M = (1.2) 

(7 R'T^y 

where 7 is the ratio of specific heats and R* the gas constant. We also define the Prandtl 
number to be 


a — 


p*c p 

~1F~' 


(1.3) 


where k* is the thermal conductivity of the fluid. 
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The non-dimensional equations of continuity, momentum and energy may then be 
written 


dp 


+ V.(pu) = 0, 


dt 

= “ V P - 777 V A A u )] + + 2 At)V-ii], 


Dt Re 

Dh Dp 1 

' m = m + $ + ^ V ' (,,VT) ' 


Re 


(1.4) 

(1.5) 

( 1 . 6 ) 


where h is the enthalpy of the fluid, and $ is the viscous dissipation. We also assume a 
perfect gas, in which case we have 


P - 


PT 

yM 2 ' 


(1.7) 


In the following sections we consider first the basic flow (section 2), which is shown to 
be a proper solution of the above equations of continuity, momentum and energy. In 
section 3 we consider the inviscid linear stability equations. In section 4 we present a 
number of numerical results, guided by which, in section 5 we develop asymptotic results 
for large wavenumbers; throughout this section we emphasise a comparison between our 
numerical and asymptotic results. In section 6 a new class of mode, which is found to 
develop at sufficiently large Mach number is considered. In section 7 we present a number 
of conclusions arising from this work. 
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2. The basic flow 

Let us consider the solution corresponding to Batchelors’ (1964) similarity solution for a 
swirling wake flow, equivalent to a far downstream limit (ar 1) of the governing equations. 
We thus suppose that the solution comprises of a uniform (freestream) flow plus a relatively 
small perturbation, i.e. 


u* 

w oo I ~ 

u = u* + 

(2.1a) 

v = 5, 

(2.1b) 

w = u;, 

(2.1c) 

U* 2 

p= ^ +p ’ 

(2.1<1) 

T = 1 + f , 

(2.1c) 

-c> 

II 

1— » 
+ 

(2.1f) 


The solution develops in much the same way as the incompressible case of Batchelor (1964) 
and it is found to leading order that 


1 



gU^logjxRe^) 

Sxu^Reoo 


Q i( r ?) + 




8 xv^Rec 


-Qiiv) 


LU^e-' \ 

Sxv^J 2 Re^ J 


( 2 . 2 ) 


5 = 0 


(2.3) 


CoRe^o f 
w = — — . , , (1 


a" 7 ? 




)■ 


where 


r 2 Re oo 


(2.4) 


(2.5) 


rr* _ e 0 2 ^ 

3 8xu* 2 Re c 


\og(xRe 00 ) + 


LU~ 

Sxu^ 2 Re c 


( 2 . 6 ) 
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and 


Qi(tj) = e T >, 


(2.7a) 


<? 2 (t?) = e 17 [log t] + ei(» 7 ) - 0.867] + 2ei(r/). 

Co and L are constants and R( x is the freestream Reynolds number defined by 


P*oo 


We also find that 


P = 


C%U, 


* 2 

oo 


8xu^ 2 U* 2 R eoo 


(1 - e ”") 2 


+ 


2ei(,)-2ei(2,)]. 


(2.7b) 


(2.8) 


(2.9) 


whilst the temperature perturbation is given by 

f =e -^(id(V)+g+ (T ~ 1)A/ ” c ° 

V 2 '4 5 ' + + 8x^ 2 B Coo 




+ 


O e ^~ 1 * eH'-V* 3 

8a — — — (- 16a- 


C 2 


C 2 


/<- 

8a7ri £^2 $)£ 2 a2iti ut 2 

"i ^ e 4 er /(£/2) — 8 — 4 — — — e 4 erf(£/2 2 )}d£ 

+ J eZ *~ J r (-2(7{-2e^ + 2 — £ 2 eV e ((£ 3 /4) 

+( 2 e^ei(( 2 / 2)) + 


( 2 . 10 ) 


where x4 and R are constants, £ — 2 7 2 , ei(z) and erf(z) are the exponential integral and 
the error function, respectively, defined in the normal manner, and M 0 0 is the freestream 
Mach number defined by 

UL 

(2.11) 


U* 

Moo = — * . , - — 

(jR'T^ U 


l J* 

= ~M. 
* 


Note that we expect M 0 0 M. We also find it useful for the remainder of the paper to 

define the lengthscale 


r = i)i = — , 

r* 

where r, = (4x/Re oa )i is the characteristic radial lengthscale. 

In the following section we consider the inviscid stability equations. 


( 2 . 12 ) 
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3. The stability equations 

We take the general basic state to be u = U(r),v = 0, w = W{r),T = To(r),p = po(r),p = 
p 0 (r) and consider small amplitude perturbations to this flow; we write 


u = U(r) + 6F(r)E + 0(6 2 ), (3.1a) 

v = iSG(r)E + 0{6 2 ), (3.1b) 

w = W(r) + SH(r)E + 0(S 2 ), (3.1c) 

p = p 0 (r)+6 7 M 2 P(r)E + O(6 2 ), (3.1d) 

T = T 0 (r) + 6r(r)E + O(6 2 ), (3.1e) 

p = p 0 (r) + Sr(r)E + O(d 2 ), (3. If ) 

where, <5 < 1, and 

E = exp[i(aar/r il + nO — actjr s ) ]. (3-lg) 


The governing stability equations, neglecting the effects of viscosity are then 

G nH\ 


T<p + p 0 ( ocF + G' + — H — —j + Gp' 0 — 0, 

po<pF + poU'G = — aP, 

„ 2 WE TW 2 

+ P0 i — P > 

r r 

Po <pH + p 0 (w' + y)c = 


/W r + doGTo(r) - 


7 M 2 P = TT 0 + rp 0 , 


7 ~ 1 
7 


( 7 <pM 2 P + Gp') = 0, 


(3.2) 

(3-3) 

(3.4) 

(3.5) 

(3.6) 

(3.7) 


where we have written 


y? = a(C7 — c) + 


nil' 


(3.8) 
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and primes denote differentiation with respect to r. Note that (3.1), (3.2) (3.7) all im- 

plicitly assume that the axial scale for the perturbation quantities is considerably shorter 
than any developmental lengthscale for the basic flow, as is the case if r 3 < Cl. 

Equations (3.2)-(3.7) may be combined to yield the following two, first order equa- 
tions 


dG 

dr 


n(rW)' + a r 2 U' p' 0 1 

7Po r 


r 2 ip 


G + 


1 ["( a 2 r 2 + n 2 ) ipM 2 ' 

tT- 


P o 


r 2 ip 


P, 


(3.9) 


dP 

dr 


= P o 




(W 2 r 2 )' W 2 ( Tq 


+ (2_1)_pL)1 

v To y 7 VoTo/J 


r 3 ip 


rip 


G 


7 ' PoTo 

[2 nW W 2 M 2 1 


L r 2 ip 


rT 0 


P. 


(3,10) 


These equations are somewhat similar to those considered by Michalke (1971), in the 
context of jet flows, if the swirl velocity is neglected. 

Let us now consider the specific basic flow of the trailing line vortex, as discussed in 

the previous section. Equations (2.1), (2.3) — (2.5) may be substituted into (3.9), (3.10), 

and then if we assume \8\ <C |Mogx|, and using the fact that a simple transformation and 

inversion of the axial velocity only affects the frequency of the stability analysis, and does 

not change the amplification factor C{ (where c = c r + ic{) and noting that p, p and T are 

an order of O(logx) smaller than u, v and then (3.9), (3.10) reduce to 

;n(rW) f + ar 2 lT _ lj r< 

r J L 


dG 

dr 


r 2 ip 


\a 2 r 2 + n 2 


r 2 (p 


-ipM 1 P, 


(3.11) 


dP_ 

dr 


(W 2 r 2 )' i \2nW W 2 M 2 ~\ 

ip ^ 

r A ip 


G - 


r 2 ip 


The boundary conditions to be applied are then 
P( 0) = 0, n + 0; P'(0), n = 0, 

G(0) = 0, |n| i 1; G'(0), |n| = 1, 


(3.12) 


and G(r), P(r) bounded as r — » oo, 


(3.13) 
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and the basic flow may be taken to be 


U = e 


W = 1(1 - e- ), 
r 


(3.14) 

(3.15) 


where q is an order one (swirl) parameter, and we have effectively scaled velocities with 
respect to U*. 

Equations (3.11), (3.12) may be combined to eliminate G , yielding the following 
second-order equation for P. 


n „ f 2 nW W 2 M 2 1 

Prr + Pr{ + - - 

I r z ip r r 

[< ar 2 U ' + n(rWy ] 


(W 2 r 2 )’ 




r°ip 


(W 2 r 2 )' 

V 3 


r 2 tp 


+ P< 


2 nW _ W 2 M 2 
r 2 p> r 


\2nW W 2 M 2 ] 

(W 2 r 2 y 1 

7 

( W 2 r 2 y 

'I 

bO 

i 

P 3 

r Y 

/ 

P 3 

r^ip 


/I _ [ar 2 U f +n(rWy] 


r 2 <p 


a 2 r 2 + n 2 


)i 


2 nW W 2 M 2 


r 2 (p 


-ipM 2 


r 2 (p 




(T^ 2 r 2 y 

r 3 (p 


= 0 . 


(3,16) 


Setting M = 0 clearly reduces (3.11), (3.12) and (3.16) to the incompressible problem, as 
considered by Lessen et al. (1974), Duck & Foster (1980), Leibovich Sz Stewartson (1983) 
and Duck (1986). It is also possible to be rather more precise regarding the behaviour of 
the solution as r — > oo; this takes the form 

dG 


dr 


±a(l-c 2 M 2 )iG = 0, 


dP r, 0,1 

— ±a l-c 2 M 2 y? = 0, 
dr 


(3.17) 

(3.18) 


where positive signs are taken for 3?{(1 - c 2 M 2 )?} > 0, and vice versa, to ensure bound- 
edness as r — > oo. 
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The “order one” problem requires a fully numerical solution, which we consider in the 
following section, prior to considering various asymptotic limits of this system of equations, 
which permit a certain amount of analytical progress. 
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4. Numerical results 

The system was treated numerically using four different techniques. The first was based on 
the method of Duck & Foster (1980), in which the system was approximated by second- 
order central differences, with conditions (3.17) and (3.18) imposed at a finite radial value 
^ ^maz, taken sufficiently large not to substantially affect the result. The determinant of 
the system was then forced to zero by adjusting the complex wavespeed c (using Newton 
iteration). 

The second method was based on a fourth-order Runge-Kutta scheme, conditions 
(3.17) and (3.18) were approximated by imposing boundary conditions at a finite radial 
value r = r max , where this value was again taken sufficiently large to not substantially 
affect the numerical results. The computations were performed by shooting the solution 
towards r = 0 and were not necessarily confined to the real r-axis. The value of c was 
adjusted (again using Newton iteration), to ensure the correct behaviour of the solution 
as r — > 0. 

The third method used was based on the first, but was a global finite difference 
method. Using the aforementioned finite— difference scheme, then by defining two addi- 
tional quantities G = c(7, P = cP at each grid point, it is possible to write the resulting 
scheme in the form 

(A — cB)x = 0, ( 4 .i) 

where 

* = [{£?}, {P), {G}, {p}] r (4.2) 

This scheme has the obvious advantage of generating 4iV eigenvalues (where N is the num- 
ber of grid points), simultaneously. The principle disadvantages are (i) it is not possible to 
use (3.17) or (3.18) because of the non-linearity of c, and consequently Dirichlet boundary 
conditions were applied instead at the outer edge of the computational domain; and (ii) 
the scheme requires rapidly increasing computational resources (both in storage and time) 
as N increases. 
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The fourth scheme implemented was a Chebyshev spectral collocation scheme, based 
on that of Khorrami, Ash h Malik (1989). This was a global method, which generally 
gave very accurate eigenvalues, but also yielded a large number of spurious eigenvalues, a 
feature found in many spectral schemes. The first (finite-difference) scheme was quick and 
robust, but because of its “local” nature, mode jumping was often experienced, due to the 
frequent close proximity of neighbouring modes, as discussed later in the paper. The second 
(Runge-Kutta) scheme was also quite fast, and had the advantage of being able to compute 
neutral and near-neutral modes (and even stable modes) by contour indentation, but 
again because it involved local searching was prone to mode jumping. The third scheme, 
namely the global finite-difference scheme, which proved to be very robust, produced few, 
if any spurious (i.e. non-physical) modes, and generated many eigenvalues simultaneously 
because of its global nature. Consequently our results were generally computed using the 
first scheme (if only a limited number of modes were required) or the third scheme if it 
was required to compute many modes. 

We now present a few numerical results to give some indication of the effects of 
variation of certain of the important parameters. Here, and indeed in all our calculations 
we chose 7 = 1.4, q = 0.8. Further, we generally found that r max = 5 was sufficient with 
A r « 0.03. 

Figures 1 a, b, c show the variation of growth rate (acj) with a for the case M = 3, 
with n = “ 1 , — 2, —3 respectively. These results (and all those presented in this paper) are 
accurate to within the graphical accuracy of the figures. Because of the great multiplicity 
and close proximity of all these modes, we show the computed values of aci at each value 
of a. From the outset it should be stated that all our results relate to negative values of 
n; in general, we believe that instabilities are almost exclusively confined to this region of 
parameter space, with perhaps a few minor exceptions. We certainly expect the largest 
growth rates to be confined to the negative values of n (this is also the case when the flow 
is incompressible). 

We see in figures 1 the following general trends: (i) an increase in growth rates as —n 
increases; (ii) an apparent cut-off value of a above which no unstable modes exist; (iii) a 
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tendency for the maximum growth rate to be attained at an increasing value of a as — n 
is increased; (iv) many modes of instability. 

In figures 2 a, b, c, d, e, growth rates for the case M = 5 with n = — 1, —2, —3, —4, —5 
respectively are presented. The trends (i) — (iv) described above are again observed, to- 
gether with the result that for corresponding n and a, the growth rates of M = 5 are 
substantially reduced compared to M — 3. Further the upper limit of a of the instability 
appears to be quite independent of M, and section 6 confirms this observation. Note that 
in figures 1 and 2, the lower range of a has been deliberately truncated. In this limit, our 
numerical results became extremely sensitive to grid size, and computation requirements 
became prohibitive. 

These results suggest a number of interesting features, and in the following sections 
we mount a systematic study for increasing M , when (— n) 1. 
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5. Asymptotic results, (-n) > 1 
5.1 M=0(1) 

It turns out that for this order of Mach number, to the orders to which we concern ourselves, 
the solution (in particular for the complex wavespeed) remains unchanged from the — n >> 1 
solution of the incompressible case, as considered by Leibovich & Stewartson (1983) and 
Duck (1986). However, since this solution forms a basis for the following subsections, we 
outline, briefly, the form of the structure in this case; full details can be found in the 
aforementioned papers. We have that 

a = na, a = 0(1), (5-1) 


and the complex wavespeed develops as 


Cl c 2 

c — Co H h 3 + • • • • 

n n 2 


Then generally 


(p2 

ip = nipo + <pi + r— T + • • • » 
n 2 


where 


<P o = a(U(r ) - c 0 ) + 


W(r) 


(5.2) 


(5.3) 


(5.4) 


However the solution is found to be concentrated about points r = ro (critical points) 
where y> 0 (r = r 0 ) = 0 and so 

a(U(r 0 ) - c 0 ) + = 0. (5.5) 

r 0 

However, it turns out that r o must in fact also be a turning point, and so 


ifj,(r 0 ) = 0, 


(5.6) 


i.e. 


cnU' [r o ) + (~) — 


(5.7) 
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(where a subscript zero here and hereafter denotes evaluation at r = ro). Equation (5.7) 
then serves to determine r 0 , and hence c 0 may be determined from (5.5). The key length- 
scale inside the critical layer is given by 


R = (r — r 0 )|n| 4 , 


(5.8) 


and then on this scale 


¥■=■¥>!- ! 7 T ^(-^2 + \<e'l{r«)R 2 } + 0(r>->), 

n 2 A 


(5.9) 


where 


= -acj, 

d 2 

<P'o( r o) = Jp-[v>o( r )] r= r 0 . 


(5.10) 

(5.11) 


For consistency it was shown by Leibovich & Stewartson (1983) and Duck (1986) that 
2 (WV)'„ [(r W)'„ + arjc/il 2 


r^a 2 


(1 + a 2 r 2 )r 2 0i 2 


(5.12) 


Further, on the R = 0(1) scale the eigenfunctions scale as 

P = P + ... , 

G = nG + ... , 


(5.13a) 

(5.13b) 


and then P(R) is described by 


Prr-{‘ 


c 2 ip'o(r 0 )R 2 1 sign(n)(l. + cy 2 rg) 


ci 


2aci 




or 





+ 


\t)C2 

2(A,)i 


}p = 0, 


(5.15) 
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(5.16) 


where we have written £ = 2 2 A * R, with 


2 

A 0 = — sign(n) 
ci 


,(1 + * 2 rg) 


and Ai = 


-Ao^o(ro) 


2d 


Equation (5.15) has a solution which may be written in terms of Weber parabolic cylinder 
functions, and so if we demand the solution decays as |£| — ■* 00 , then m must be 

an integer, yielding the following result for c 2 : 


C2 = 


(A.)’ 

Ao 


(1 + 2m), 


m = 0, 1,2, . . . 


(5.17) 


Figure 3 shows the variation of growth rates for M = 3, for n — — 1 (the least unstable 
mode shown) up to n = — 10 (the most unstable mode shown) as computed for the full 
system; in all cases the most unstable mode for each value of n is shown. Using the 
asymptotic results above, we show the corresponding results in figure 4 (in particular we 
set m = 0 in (5-17)). The comparison between figures 3 and 4 reveals good agreement in 
the growth rates at larger values of — n, although at smaller values of — n, comparison is 
less good than the comparison of incompressible numerical and asymptotic results. This 
less satisfactory agreement may be attributed directly to the effects of compressibility, but 
for |n| sufficiently large, the asymptotic results presented in this sub-section are ultimately 
approached. 

In the following sub-section we begin to incorporate compressibility into our asymp- 
totic description. 


5.2 M=0(|n|l) 

The key equation (5.14) was obtained by taking the 0(|n|«) terms in equation (3.16). 
If |n| 1, and we now permit M to grow in magnitude, additional terms will ultimately 

enter (5.14) when 

M = M = 0(1), ( 5 - 18 ) 
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when the coefficient of P will include the effects of compressibility. The modified equation 
is then 


Prr ~ sign(n) |2 [— - 
t Lei 


2cvcj 

+ 


2 

(1 +« 2 r 0 2 ) 


■ 

r o 


r (rW); + arjF{ 1 

(It'.V)! 


occirl J 

\ r 0 / / 


)}p = 0. 


(5.19) 


All other quantities (specifically cq and ci) remain unchanged from those evaluated previ- 
ously. This equation may be transformed to the same form as (5.15), viz 


« 


_{l f+ A^ i _A 2 )w = 
14 2(Ai)* J 


(5.20) 


where Ao and Aj are given by (5.16), whilst 


A 2 = sign(n) 


(rwy 0 + arjUp 

drgci 


WqM 2 


r 0 


(5.21) 


Consequently, using our previous arguments 
(^ 1 ) 2 (1 + 2m) + A 2 


C2 




Ao 




m — 0, 1, 2, . . . 


(5.22)) 


Figure 5 shows a comparison of “exact” growth rates (obtained using the numerical ap- 
proach of section 4), shown as a broken line, with the asymptotic results of this subsection; 
here, the chosen values of n are -5, -10 and -15 and we set M = |n|? i.e. M = 1 in each 
case. The agreement is satisfactory, and improves as n increases. Figure 6 shows the 
corresponding comparison for c r ; in this case the agreement is excellent. 


5.3 M=0(|n|r) 

As M (and hence M) increase in magnitude, the “A 2 component” of C 2 increases (as 
M 2 ) and becomes sufficiently larger than “Ai component”. Simultaneously, the coefficient 
of P r in (3.16) will grow, and further terms in the coefficient of P will become significant. 
The next important regime for M is when 

M = |n|*M, M = 0(1). (5.23) 
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Although the key radial scale remains R = 0(1) (see (5.8)), the series development of c 
and ip(R) is now altered, and is instead 


c i 

c — Co H — h 

n 


C 2 


Ini' 


C3 

|n|f 


+ ... 


(5.24) 


and therefore, 


¥ ’ = v ’ , + nT + wT + - 


= —aci - -~-ac 2 sign(n) 
n 7 


-^|-ac 3 sign(n) + ^o(r 0 )i? 2 } + ... . 


(5.25) 


Considering 0(|n|<) quantities in (3.16) requires that the sum of these terms is zero, and 
so this leads to, 


2ar 0 


C2 = -si g n(n) M2 W »r (r ^ + ^ 


1 + a 2 r 2 


(5.26) 


The equation for P(R ) is again obtained by taking the 0(|n|2) terms in (3.16), namely 


Prr 


MHV b R - { 1 + f 1 i f 2 ^i + } }p = o. 

I ri L ip i J J 


ro 


(5.27) 


Using the standard transformation 


P(R) = P{R) exp 


IM^WqR] 

-2 r 0 


then 


- R -{ 2signW (|-«)( li ^) 


s c\ 2ac\ 




(5.28) 


(5.29) 
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or 


where 


U 2 + — C2+ , A3 P 


2 (Aj ) 3 


P = 0 , 


(5.30) 


M 2 Wq 


(5.31) 


Since (5.30) is again a form of Weber’s equation, then 


~{ 


(Ai)^(l + 2m) + A3 


}, m = 0,l,2 ,... 


(5.32) 


Figure 7 shows a comparison of asymptotic results (solid lines) with numerical results 
(broken line), for n =-5, - 10 , -15 with M — |n|« (i.e. M = 1 ). Again the comparison 
becomes markedly better as |n| increases. The corresponding distribution of e r is again 
good, and shown in figure 8 . 

5.4 M=0(|n|i) 

In this case it is quite clear from the previous subsection that when M = 0(|n|s), then the 
C 2 term (see (5.26)), which grows as M 2 will become comparable in magnitude to the term 
involving ci (which is independent of M ); in addition to this, since C 3 grows as Af 4 as M 
increases, then this term will also become comparable to the c\ term when M = 0(|n |?). 
In this case, it is straightforward to show that R = 0(1) remains the appropriate radial 
scale, whilst the wavespeed expands as follows 


_ , Cl , C2 
c — Co H b - — -3 + • • • 

n n 2 


(5.33) 


implying 


¥> = <Pi + - V + -.. • 

h * 


(5.34) 


If we write 


M = \n\*M, M = 0(1), 


(5.35) 
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r o 


then the equation for P(R) is 

l LV rfoi Ar jjv?i r n / 

'+^o\/ (^ 2 r 2 )[, \-| 

J’oV’i 

'(rW)' Q +arlU',\^ 2 


+ 


+ 


- 2 (> 


*o¥>i 


H> i 


+sign(n) ( 


(rb^+Ar 2 ^ ^ M 2 ^ ^ 




+2 


To ' V 3 ! 


/l + a 2 r 2 w(^ 2 r 2 )' 


V r 2 v?! ) V r 2 v?j 

) ipx\) 


Again using a standard transformation of the form 

’1 \n\< M 2 Wq R) 


P(R) = P(R ) exp 


ro 


leads to the differential equation 

liU (rWY 0 +arlU^ f 2W a 


P RR — jM 2 ( 


+ 


r lVl 


\ rfa i / V 

(rW0;+ar=E7JVv; 2 


0 ( 




+ 


>( 


-sign(n)— ) 

ro ' 

(trV) ^ , [tiiwsy] 

r aV\ 


0 + (^) 1 


r lvl 


-y 


¥>1 


' rfo i /V r 0 / 


/l + d 2 r 2 \ 

{ (W 2 r 2 )' 0 


V r;fci ) 

' ^¥>1 

' V 7 ! d 


(5.36) 


(5.37) 


(5.38) 


If our assumption regarding the importance of the R = 0(1) scale is to be consistent, 
(although see the comments below regarding R = 0(|n|4)) then the coefficient involving 
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|n| 2 above, must be zero. (This is also consistent with the previous smaller orders of Mach 
number considered previously). This leads to 


, (rW)' 0 +arlVl } // A/*H'? y (W^r*)i,r 0 ~ 

2dr 0 \ 1 + <x 2 r1 ) y \ 2 oTq / l + o: 2 ?-^ C * 

— - 

+ 1 + « 2 r2 V 2 / 


(5.39) 


where Cj nc denotes the incompressible value of Ci, given by (5.12). The equation that 
determines P(R ) is then, 




. , x WqM ( 

+s!gn ( „ ) __( 


Wq M 2 f(rW)'„+„rlV;, 


rlac i 


)ir 


0, 


(5.40) 


or symbolically 

P RR — + ^ 0 ^ 2 ] P — 0, 


(5.41) 


and using previous arguments, we must have 

i)*(l + 

Ac 

Note that although the transformation (5.37) suggests growth a s R —> 00 , this is more than 
offset by the decay of the parabolic cylinder functions, albeit on a larger lengthscale R = 
0(|n|4). Note too that the transformation (5.37) is consistent with that used previously, 
namely (5.28). 

Figure 9 shows comparison between the above asymptotic results (solid line) and the 
“exact” numerical results (broken line), for the most unstable growth rates when n = -5, 
-10, -15 (M = |n|^, i.e. M = 1 in all cases). Again the agreement is seen to improve as 
— n increases. Figure 10 shows a comparison between the corresponding c r ’s, for the above 
cases and indicates good agreement between our asymptotic and numerical results. 


2m) 


m = 0, 1, 2 . 


(5.42) 
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It turns out, however, that this order of M marks a watershed. When M is not 
large, both roots of c\ must be complex (and from our previous discussions, the flow is 
unstable). As M increases, however, eventually terms inside the square root term in (5.39) 
will eventually become positive, and hence the two roots of c\ will cease to be complex 
conjugate pairs, but will both become real. In particular, this will occur when 


MjW 0 2 (W 2 r 2 )'r 0 

2 ar 2 1 + a 2 TQ 


(5.43) 


To illustrate the stabilisation of these modes, in figures 11 a, b, c, we show the variation 
of the growth rates with M for the cases n = —5 and a = 2.5, n = — 10 and a = 5, 
n = — 15 and a = 7.5 respectively. We clearly see these modes becoming neutrally stable 
at finite values of M. For comparison, in figures 12a — 12c (corresponding to figures 11a — 
11c respectively) we show the corresponding results from solutions of the full system (3.11), 
(3.12). At the lower values of M, there is good correlation between the two sets of results. 
However, as the Mach number increases, and as the order of the modes (i.e. m) increases, 
the correlation deteriorates. We consider this latter point first by examining the behaviour 
of the higher order modes. 


5.5 The structure when m =0(|n|*) , (M=0(|n|^)) 

One detail that has not been considered so far, is the behaviour of the modes as m increases. 
This feature was considered by Duck (1986), and we consider this aspect here, for the 
important case when M = 0(|n|^). c develops in the same manner as in the previous 
sub-section, and close to the critical layer y takes the form 

V 3 - VI + TTTV^ + • • (5.44) 

|n| 2 

co remains unchanged from that given previously, however Ci will differ. In this case, we 
consider the lengthscale 


R = (r - r 0 )|n |2 = 0(1), 


(5.45) 
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and so on this scale 


Wi = —OiC\ + 
If 2 = — «C 2 + 


sign(n) 


3! 


<Po(r 0 )R 3 . 


(5.46) 

(5.47) 


Substitution of these expansions into our governing equation gives, to leading order, an 
eigenfunction equation of the form 


P M + q,(R)P^+\n\q 2 (R)P = 0- 


R 


(5.48) 


Using the standard transformation 


P(R) = P*(R)exp 


i(*P 


(5.49) 


we obtain 


P-Z + [\n\q 2 ~ R~ A) P * -°- 


(5.50) 


However, since |n| is large, (5.50) may be approximated by 


Pi g +|n| ?2 P*=0, 


(5.51) 


where 


92 = 


r oV?i 


4r„Wg + (1 +a 2 r 0 2 )(rj v ; - 


-sign(n)2Wo r 0 M <p 1 + 


3 1 


(5.52) 


Thus (5.50) may be written as 


p; l + h > 




[ (l + t 2 ) 0 + * 2 ) 


P * = 0 , 


(5.53) 
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with 


2[4r 0 W 0 2 - (1 + a 2 rg)(iy 2 r 2 )o]sign(7i) 
® c i r oV ? o( r °) 

— 2[(1 + d 2 r 2 ) + MHg/£acisign(n) 
r o¥>o( r o) 

AW*M 2 

r o^o( r o)’ 


£ = 


-y , o(?~o)sign(n) 

2ac\ 


1 

2 




(5.54) 

(5.55) 

(5.56) 


(5.57) 


h 2 = |n|. 


(5.58) 


Equation (5.53) is now in a form suitable for a WKB J type of approximation. This equation 
has four turning points, however for the range of £ required that we need consider only 
the turning points at 

1 

2 

(5.59) 


£ = ± 



+ -^-1 


For large h the WKBJ solution is given by 


P* = q 


A\ exp J (ihq? d£) + Ai exp J —(ihq' 2 d£) 


(5.60) 


where 


(1 + £ 2 ) 2 0 - + £ 2 ) 


(5.61) 


The treatment near the turning points is standard (see for example Duck 1986) and leads 
to the following dispersion relationship for c\ 


Ih = — — (1 + 2m), m = 0,l,2,... 


(5.62) 
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where 



(5.63) 


6 = 


1 77 1 

+ ~ 

2u 2 


1 



(5.64) 


9 = ~ 


(S + 5f) 1/4 ’ 


9 1 2 ’ 

r 


k = g [ — — + — <7 1 — 1 ) , 

' 2u 2 ' 


(5.65) 

(5.66) 

(5.67) 




- 2fc 2 

-9i)g 2 ' 


(5.68) 


Here, ^(y,^), E(^,k), n(^,a 2 ,Ar) denote complete elliptic integrals of the first, second 
and third kinds respectively. 

The system (5.62) was solved using Newton iteration, and results for variations of 
growth rate with M were computed for m = 0,1, 2, 3, 4, 5 for the cases n = —5, a = 2.5 
(figure 13a), n — 10, a = 5 (figure 13b), n = —15, a = 7.5 (figure 13c). The modified 
theory of this subsection does indicate some important improvement in the comparison 
with figures 12; in particular the “bunching up” of the modes with an increase in order is 
captured. However, although for n — —15 there is quite good correlation for the lower order 
modes, as m increases, the various wiggles observed in figures 12 are not described, and 
more importantly our asymptotic results do not capture the instability shown in figures 12 
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beyond about M = 8. Indeed, figure 12 indicates instability as M increases. We therefore 
conclude that the nature of these modes, as M — > oo is somewhat different. It turns out 
from the analysis of the following section that as M — » oo, a further class of mode emerges, 
quite distinct from those considered so far. 
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6. Centre Modes, M=0(|n|) 

Here we examine modes exhibited when the Mach number is of order n, i.e. 

M = \n\M, M = 0(1). (6.1) 

As noted previously, as M —> oo, the features of the modes exhibited at lower Mach 
numbers cease to exist. In fact the analysis of sub-section 5.5 suggests that the modes 
in this regime under present consideration are neutral. However, our numerics belie this 
and point towards the existence of centre modes, by which we mean, the eigenvalues are 
determined primarily by conditions close to the axis of the vortex, r — 0, in a manner 
similar to that of Stewartson and Brown (1984, 1985). Here, the complex wavespeed 
develops as 

c — c o H — 2 * H“ • ■ • ’ (6-2) 

n z 

(an expansion that can be verified a j posteriori) and therefore 


< P ~ nip o H b . . - 

n 


= n{oc(U - c 0 ) + + {— aci + ^v?o(0) r2 } 


+ ... 


In order to be consistent we must have that tpo(r = 0) = 0, i.e. 


co = 1 + —• (6-4) 

a 

We must now go on to find c\ in our complex wavespeed expansion. The solutions are, 
in the main, trapped in the region around r = 0, but we begin by showing how these 
solutions are connected to the outer flow. The flow is divided into the four regions, which 
are considered in turn. 

We first begin our analysis at the outer region of the flow, where r = 0(1), and the 
governing equation has the form 


W 2 M 


P r + \n\ 4 M (plP = 0. 


Using a transformation of the form 


P = 



n 2 W 2 M 2 , 1 

— ' lr j' 


and writing 


Xo = -M 2 ^ o + 


W 4 M 
4r 2 ’ 


reduces (6.5) to 


P’rr - |n| 4 x „P* = 0. 


(6.6) 


(6.7) 


( 6 . 8 ) 


We define the critical point fo by xo(fo) = 0 and when r > fo, Xo < 0 and the approximate 
solution to (6.8) is given by 


P* - 


E , 


•Xo]* 


exp 


{?|n | 2 f [-Xo]’drj, 

J Tq 


(6.9) 


where E\ is an unknown constant. 

We have also assumed that [ — Xo] ^ is slightly imaginary, with its imaginary part 
positive. This condition is to ensure that the solution is bounded as r — > oo, otherwise if 
/m{[— xo]^ } < 0 we replace T by ‘—V in (6.9). Looking at this point from a more physical 
perspective, we require waves to propagate out from the critical layer and not in toward 
it. 

In order to be able to match (6.9) with the solution in the r < fo region, we must first 
examine the solution of the flow in the transition layer, around r = fo. The lengthscale in 
this region turns out to be 


r — f o = -R|n| » , R = 0(1), 


( 6 . 10 ) 


and use of this scaling reduces our governing equation to the Airy equation, which has 
solutions of the form 

P' =F,Ai(x;(r'o))ifl) + P 2 Bi( x ;(7o))*fl). (6.11) 
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Standard analysis and matching with (6.9) yields the result 


Fi = iFo. 


In the region r < fo, Xo >0 and hence the WKBJ solution here is given by 


P* = —\D i exp 

Xo 


a: 


r W 4 M -i-j -2 


n 


4r 2 


-M 2 ip 2 0 5 dr | 


H 


+D 2 exp<j - / mi 

*ro 


rW 4 M TjUjii. .. 


As we approach the transition-layer this solution must match with (6.11) implying 
Z?i = — iD 2 V3. 


Note that as r — * 0 

P* ~ Jr T <{ D,Zi exp( 

/7 A yfr 2 I \ 

where, 


^J B|Ziexp (i!^£^W!) + D 2 Z 2 exp(. r2|n|2 

2 ■ \ 4 / V 


). 


Z\ = exp{ f |n| 2 (xo 

J Tq 


q 2 M r 


) dr ■ 




We now consider the lengthscale, r = 0(|n| 2 ), by setting 
f = r|n| 2 = 0(1). 

The governing equation in this region, to leading order, is given by 


— \n\q l M rPf + |n 


1 |2 

r2 q 3 M 2 1 I 

|n| 

■ <Pi r 2 . 


P = 0, 


( 6 . 12 ) 

(6.13) 

(6.14) 

(6.15) 

(6.16) 

(6.17) 

(6.18) 
(6.19) 
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where 


¥>i = ^o(°) r2 > 


( 6 . 20 ) 


and as previously, it is necessary to use a transformation, which is of the standard type, 
i.e. 


P = P** exp 


{| J \n\q 2 rMdr |, 


which leads to the final form of the governing equation in this region, viz 

r2 


P rV + M 


2 q 3 M 1 q*r 2 M 
~p2 


^o(0)r 2 r 


4 


P** = 0 . 


This equations yields solutions of the form 

p** = r-i ( Ct ) + cyr.( ] , 


Now, as r — > oo 




\J~2TUp- 


exp 


M |n|r 2 /4 


(q 2 M |n|r 2 ^ 


( 6 . 21 ) 


(6.22) 


(6.23) 

where I and K are modified Bessel functions, written in standard form, with 

i/ = M(i (6.24) 

2 v 


/ * rml 

f g 2 M 2 |n 

If! 1 ; 

j 2 q 2 M |n|r 2 /4 

k 4 ; 


), (6-25) 


and comparing this with equation (6.15), we see that in order for these solutions to match 
and taking into account equation (6.14), then 


C x = -i I ttV3Z^C 2 . 


(6.26) 
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Also, we take note of the fact that as r — ► 0, then 


P** 




1 _ 

r 2 


v r(*/ + i) 



(C.27) 


2 

where i = q 2 M |n|/4 and r(i/) is the gamma function, defined in the standard fashion. 
Next we set 


R = r\n\ = 0(1), 


(6.28) 


and in this regime ip develops as 


¥ 


<Pi 


n 


+ ... 


= K“ 5ci + 5^'(° )JJ2 ) 


+ 


(6.29) 


with the governing equation 

4a 2 ^ 2 + 4 aq 


Prr + M 


which has the solution 


+ 


2q 3 M 


ac, + \p>q(Q)R 2 } 2 [-ac, + ^'(O)^ 2 ] 


1 ' 

^R 2 . 


P = 0, 


(6.30) 


P = 


(T R 2 )^ + < 


(1 - T R 2 ) 


j-F(wi]xv 2 ’, M + 1; Tir), 


(6.31) 


where, F is the hypergeometric function in standard notation, and 

T M(Q) 

ac\ ’ 

= ^(1^1 — *5) — /i, 

^2 = ^(M + s) - A, 

s = (\n\ 2 + b)K 

A(A + 1) = 


(6.32a) 

(6.32b) 

(6.32c) 

(6.32d) 

(6.32e) 
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with 


2 |4so± 4^\ (6.33a) 

1 Krt'(o)) 1 

V|n| 2 M 2 (6.33b) 

rt'(o) ' 

Note that this equation is very similar in form to that found by Stewartson & Capell 
(1985). This solution is of course finite at R = 0 by the definition of the hypergeometric 
function, and therefore satisfies our boundary conditions at the centre of the vortex. As 
R becomes large 

P ~ R^AxR 3 + A 2 R~ s ), ( 6 - 34 ) 


matching this with (6.27) and making use of (6.26) we see that, 


T(s/2 + l)r(s/2)H» 

(iA)’Z? 


— inA22y/3i 


(6.35) 


which gives us a relationship between A\ and A 2 > We now go on to determine c\ explicitly. 
We begin, by writing P, as R — * oo, in terms of its asymptotic expansion in this 

limit, i.e. 


(TP 2 ) 1 


(1 - T R 2 Y 


<o(-t R 2 r 


.M 

2 


+■ 




-T R‘Y 


2 


S E, 


(6.36) 


where, 


Co 



|n|! (s - 1)! 



,):(■ 


n\ s 

— H h u 

2 2 H 


y 


So 


n|! (— s — 1)! 



(6.37a) 


(6.37b) 
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and Si an d Yj 2 are k Qt h power series expansions in 1/T R 2 , with their leading order terms 
unity. Matching the two large R solutions, (6.34) and (6.36) yields the relationship, 


Ai 

A 2 


S(-r y. 

Oq 


(6.38) 


However, since |<s| >> 0 then |T 5 | 0, which in turn implies that |co I |£o|, and therefore 

to leading order, <^o = 0. Hence, we see that, ^(|n| -f s) — fi must be either a negative 
integer or zero, i.e. 


fio = g(l n l + s ) + N* 


(N > 0 ), 


(6.39) 


for integer N, where /to denotes a first order approximation to ft. 

Since s is imaginary, the above expression substituted into (6.32e) and (6.32a) gives 
the leading order imaginary contribution to cj . i.e. it is obtained from 


q^o(Q) 

8aft 0 (p-o + 1 ) 


(6.40) 


It is also possible to obtain the correction to this term by retaining the next order approx- 
imation to fi. On setting, fi = fto -f A fi, where [A/r| <C 1, then we obtain 

(-l) N+1 |n|!iV!(s - 1)! 


Co 


( s -f |n| + A)! 


-A fi, 


(6.41a) 


6 0 = 


(~l) N \n\l(s + N)\ 
s\ (|n| -f- N)l 


(6.41b) 


and then (6.38) gives 


Art = — ( + l w l + N)l(s + N)\ 

' A 2 K J s\ (s — 1)! N\ (|n| + N)\ 


(6.42) 


Since 


fi{fi + 1 ) = 


^o(O) 

4acj ’ 


(6.43) 
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and therefore, if we allow jj, = fio + A ft, in the above, we see that 


C] = 


av?o(0) 


1 - . (2/i . + 1) -. A I . 


8ajj,o(fio + 1) \ 4/io(Ao + 1) 


Therefore, we may now write c\ as 

a V 7 o (0) 


c 1 


8afi 0 (fi 0 + 1 ) 


1 + 


Aj 


(s + \n\+N)\(s+N)\ 

TA 1)1 s!(s-l)!iV!(|n| + iV)! 


X 


(2//q + 1) 

4/io(Ao + 1) 


(6.44) 


(6.45) 


where, all variables in this expression have been previously defined. However, the correction 
term is found to be of several orders of magnitude smaller than the leading order term of 
ci- 

In figure 14a we show the variation of the growth rate with a for the case n = — 5, M = 
5. This figure is to be compared with figure 2e, computed from the full system; comparison 
between the asymptotic and fully numerical results is surprisingly good, considering the 
smallness of -n and M. Figures 14b shows the growth rate variation, predicted from the 
analysis of this subsection, for the case n = — 10? Af = 10 and is to be compared with the 
fully numerical results shown in figure 15. Again, the results are encouraging, in particular 
the magnitude and location of the maximum growth rate is quite accurately predicted. 
Unfortunately, our attempts at a numerical solution of (3.11), (3.12) with a reasonable 
number of grid points for the case n = -15, M = 15 proved unsuccessful, due to a large 
number of spurious (i.e. highly grid-dependent) modes. Indeed, it was found that as M 
was increased, our numerical scheme often failed abruptly, with a sudden generation of 
large numbers of spurious modes. 

Referring again to figures 14, we see that the upper neutral point is clearly seen. This 
occurs at a = — n/q , the point at which a = 0 (see (6.33)). The fully numerical results 
throughout this paper all very clearly show this result, which appears to apply universally, 
as in the incompressible case. 
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Additionally figures 14 show a distinct lower neutral point. This corresponds to the 
value of a where <^q( 0) = 0 namely a ~ —qj 2. At this location, a and b (defined by 
(6.33)) both exhibit a singularity. There is some correlation with these results and those 
of the fully numerical scheme. Indeed, it is most remarkable that the location of this 
lower neutral point is identical with that found in the incompressible case (Stewartson & 
Leibovich 1987), and that both the upper and lower neutral points are captured by (6.45). 

It is also worth noting that there is some similarity between the structure of these 
modes, and those found in the incompressible work of Stewartson & Brown (1985), al- 
though the particular details are different, and in our case analytic/asymptotic solution on 
the r — 0(1) scale is possible (see 6.6), whilst Stewartson & Brown (1985) had to resort 
to a numerical approach for this scale* 

In the following subsection we go on to draw a number of conclusions from our study. 
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7. Conclusions 

We have mounted a systematic study of the inviscid stability of the trailing line 
vortex, starting at zero Mach number M, and progressively increasing M. We see a general 
reduction in growth rates as M increases, and indeed the results of subsections 5.4 and 
5.5 predict that the original family of modes will stabilise when M = 0(|np), specifically 
when (5.43) is satisfied. However it is shown in section 6 that when M = 0(|n|), a centre 
mode class of instability is formed. We feel that although our numerical results (figures 
12 in particular) indicate that these modes spring from the higher order modes at lower 
Mach numbers, as — n increases these centre modes modes may well become distinct from 
the original class of modes that exist at lower Mach numbers. 

A further important feature of note, and one that is observed in incompressible work 
(Leibovich Sz Stewartson 1983, Stewartson & Brown 1985) is that there is a good deal of 
numerical evidence to suggest that there exist no instabilities for a > —n/q. Indeed, this 
is also confirmed by (6.45), noting that a is given by (6.33). This aspect is currently under 
further investigation. 


35 



Acknowledgement s 

J.A.K.S was in receipt of a S.E.R.C. studentship. A number of the computations were 
performed with the computer facilities provided by the University of Manchester. 


36 


References 

Batchelor, G.K. 1964 Axial flow in trailing line vortices. J. Fluid Mech. 20, 645 

Coleman, C.S. 1989 The stability of swirling jets. Astro. Soc. Australia Proc. 8, 38 

Duck, P.W. 1986 The inviscid stability of a swirling flow: large wavenumber disturbance. 
J. Appl. Math. Phys. (ZAMP) 37, 340 

Duck, P.W. & Foster, M.R. 1980 The inviscid stability of a trailing line vortex. J. Appl. 
Math. Phys. (ZAMP) 31, 524 

Duck, P.W. & Khorrami, M.R. 1991 On the effects of viscosity on the stability of a 
trailing line vortex. ICASE report 91-6 (To appear in J. Fluid Mech.) 

Foster, M.R. & Duck, P.W. 1982 Inviscid stability of Long’s vortex. Phys. Fluids 25, 
1715 

Foster, M.R. & Smith, F.T. 1989 Stability of Long’s vortex at large flow force. J. Fluid 
Mech. 206, 405 

Gortler, H. 1954 Theoretical investigations of laminar boundary-layer problems II — the- 
ory of swirl in an axially symmetric jet, far from the orifice. Air Force Contract no. AF 
61(514) 627C 

Khorrami, M.R. 1991 On the viscous modes of instability of a trailing line vortex. J. Fluid 
Mech. 225, 197 

Khorrami, M.R., Malik, M.R. & Ash, R.I. 1989 Application of spectral collocation tech- 
niques to the stability of swirling flows. J. Comp. Phys. 81, 206 

Leibovich, S. & Stewartson, K. 1983 A sufficient condition for the instability of columnar 
vortices. J. Fluid Mech. 126 . 355 

Lessen, M. & Paillet, F. 1974 The stability of a trailing line vortex. Part 2. Viscous 
theory. J. Fluid Mech. 65, 769 


37 



Lessen, M., Singh, P. & Paillet, F. 1974 The stability of a trailing line vortex. Part 1. 
Inviscid theory. J. Fluid Mech. 53, 753 

Loitsyanskii, L.G. 1953 Propogation of a rotating jet in an infinite space surrounded by 
the same liquid. Prik. Mat. Mehk.17, 3 

Michalke, A. 1971 Instabilitateines Kompressiblen rundens Freistrahls inter Beruck- 
sichtingung des Einflusses den Strahlgrenzschtdicke. J. Flugwiss 19, 319 

Michalke, A. 1984 Survey on jet instability theory. Prog, in Aeropace Sci. 21, 519 

Stewartson, K. & Brown, S.N. 1984 Inviscid centre-modes and wall-modes in the stability 
theory of swirling Poiseuille flow. IMA J. Appl. Math. 32 , 311 

Stewartson, K. & Brown, S.N. 1985 Near neutral centre modes as inviscid perturbations 
to a trailing line vortex. J. Fluid Mech. 156 , 387 

Stewartson, K. & Capell, S. 1985 On the stability of ring modes in a trailing line vortex: 
the upper neutral points. J. Fluid Mech. 156 . 369 

Stewartson, K. & Leibovich, S. 1987 On the stability of a columnar vortex to disturbances 
with large azimuthal wavenumber: the lower neutral points. J. Fluid Mech. 178 , 549 


0.06 



Figure la Variation of growth rate with a, n = —1, M = 3. 



Figure lb Variation of growth rate with a, n = —2, M = 3. 
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Figure 2a Vi 
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Figure 11c Variation of growth rate with M using section 5.4 results, 
first 6 modes. 



Figure 12a Variation of growth rate with M fully numerical results, n 
first 6 modes. 
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Figure 12b Variation of growth rate with M fully numerical results, n = —10, a 
6 modes. 
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Figure 12c Variation of growth rate with M fully numerical results, n 
first 6 modes. 
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Figure 14b Variation of growth rate with a for centre modes, n = —10, Af = 10. 
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Figure 15 Variation of growth rate with a n 


— 10, M = 10, fully numerical results. 
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